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ABSTRACT 

Neutron stars spinning at moderate rates (^300-600Hz) become oblate in shape and acquire a 
nonzero quadrupole moment. In this paper, we calculate profiles of atomic features from such neutron 
stars using a ray-tracing algorithm in the Hartle-Thorne approximation. We show that line profiles 
acquire cores that are much narrower than the widths expected from pure Doppler effects for a large 
range of observer inclinations. As a result, the effects of both the oblateness and the quadrupole 
moments of neutron stars need to be taken into account when aiming to measure neutron star radii 
from rotationally broadened lines. Moreover, the presence of these narrow cores substantially increases 
the likelihood of detecting atomic lines from rapidly spinning neutron stars. 
Subject headings: stars: neutron — relativistic processes — gravitation 



1. INTRODUCTION 

Atomic features from neutron star surfaces, such as ab- 
sorption lines, emission lines, and edges, provide some of 
the most direct probes of the physical properties of these 
objects. The origin of the features depends on the com- 
position of the stellar atmosphere and the conditions at 
the stellar surface, such as the temperature, the metal- 
licity, and the magnetic field strength. 

Apart from probing the conditions of the stellar sur- 
face, it has long been recognized that atomic features can 
help constrain the neutron star radius R and its mass M. 
In particular, two global parameters of a neutron star can 
be measured from observations of such spectral features. 
The redshift of the line center can be used to determine 
the compactness 2GM/Rc 2 of the neutron star. At the 
same time, the line width for a neutron star of known 
spin frequency can lead to a measurement of its radius 
(Ozel & Psaltis 2003). Together, these two parameters 
can constrain the equation of state of neutron star inte- 
riors. 

The detailed profile and overall redshift of atomic fea- 
tures from neutron-star surfaces is dictated primarily by 
the properties of their spacetimes. The spacetimes of 
slowly spinning neutron stars can be approximated by 
the Schwarzschild metric, with the addition of Doppler 
shifts due to the velocity of their surfaces (e.g. Miller & 
Lamb 1998; Muno et al. 2002; Poutanen & Beloborodov 
2006). In this approximation, the rotational broadening 
of spectral features is largely due to the Doppler shift and 
is approximately proportional to the neutron star radius 
(Ozel & Psaltis 2003; Chang et al. 2006). 

The Schwarzschild approximation neglects three effects 
which become important for more rapid rotation. To first 
order in stellar spin, the rotation of the star causes frame 
dragging in the spacetime around it. Bhattacharyya 
et al. (2006) showed that the effects of frame dragging 
on spectral lines are small for the 300-700 Hz spin fre- 
quencies observed from X-ray bursters and millisecond 
pulsars. However, stars spinning above 300 Hz become 
oblate in shape and acquire a non-zero quadrupole mo- 
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ment (e.g., Cook et al. 1994; see also Laarakkers & 
Poisson 1999; Morsink et al. 2007). For such stars, 
the Schwarzschild+Doppler method is no longer accu- 
rate (Chang et al. 2006). Instead, a different metric is 
needed that takes into account higher order effects on 
the neutron star spacetime. 

For moderately spinning stars (< 800 Hz), the exter- 
nal spacetime can be approximated analytically to sec- 
ond order in spin frequency by the Hartle-Thorne metric, 
which allows for an arbitrary quadrupole moment (Har- 
tle & Thorne 1968). The spacetimes around more rapidly 
rotating neutron stars can only be modeled with numer- 
ical metrics (see, e.g., Cook et al. 1994; Stergioulas & 
Friedman 1995). 

In this paper, we calculate the profiles of atomic fea- 
tures from neutron-star surfaces in the Hartle-Thorne ap- 
proximation. This approach allows us to describe the ex- 
ternal spacetimes of rotating neutron stars using only a 
few of their macroscopic properties, such as their masses, 
radii, angular momenta, and quadrupole moments, with- 
out a detailed knowledge of the equation of state; nat- 
urally, these properties can be calculated for any given 
equation of state. We find that taking into account the 
quadrupole moment of the neutron-star spacetime leads 
to atomic features that are qualitatively different than 
those calculated in the Schwarzschild+Doppler approxi- 
mation. Indeed, atomic features with very narrow cores 
may originate on the surfaces of even moderately spin- 
ning neutron stars when viewed at inclinations at which 
pure Doppler effects would predict significantly broader 
features. 

2. RAY TRACING 

We describe the external spacetime around a neutron 
star by a variant of the Hartle-Thorne metric (Hartle 
& Thorne 1968) constructed by Glampedakis & Babak 
(2006). We write this metric in Boyer-Lindquist coordi- 
nates as a deviation from the Kerr metric: 

= 9^ + Va 2 h^, (1) 
where g~ v is characterized by the line element 
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Here M is the mass of the neutron star and a its specific 
angular momentum per unit mass. 

The parameter rj measures the deviation of the 
quadrupole moment, q, of the neutron star from the Kerr 
value 



q = -a z (l + r 1 ). 



(5) 



In the Kerr metric, rj — and Equation ^ reduces to 
q = —a 2 . The correction to the Kerr metric is given by 

h tt = {\-2M/r)- 1 [(1-3 cos 2 6») Ji(r)] , 
h rr = (1 - 2M/r) [(1 - 3 cos 2 6) 7i(r)] , 
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where the functions J-1,2 are given explicitly in Ap- 
pendix A of Glampedakis & Babak (2006). 

Laarakkers and Poisson (1999) and Pappas & Apos- 
tolatos (2012) calculated the value of r\ numerically for 
several equations of state and spin frequencies. They 
found that Equation ([5| provides a good description of 
the quadrupole moment even for stars spinning near their 
breakup frequency and that, in general, the parameter 
r] varies between 1 and 6. For the purposes of the il- 
lustrative calculations presented in this paper, we have 
adopted r\ — 3.3, roughly corresponding to a 1.4 M star 
with an FPS equation of state (Laarakkers & Poisson 
1999). We have also adopted a value of a = 0.357 for 
a spin frequency of 700 Hz. We adjusted this value lin- 
early with spin frequency in the calculations for slower 
neutron stars. 

Neutron stars spinning at moderate to fast spin fre- 
quencies also acquire an oblate shape (Cook et al. 1994). 
In our calculations, we allow for the surface of the neu- 
tron star to deviate from spherical in addition to adopt- 
ing a non-zero quadrupole moment. Morsink et al. 
(2007) approximate the oblate shape of the stellar sur- 
face by expanding the radius R{6) as a series of Legendrc 
polynomials: 
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where Pi n is the Legendre polynomial of order 2n and the 
coefficients ai n can be represented by simple polynomial 
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Fig. 1. — Line profiles for neutron stars spinning at three different 
rotational frequencies. In each case, the star has a radius of 10 km, 
a mass of 1.4 Mq, and an inclinations of (top) 90° and (bottom) 
20° to the observer. As the spin frequency increases, the lines 
become broader. At low observer inclinations, the lines also develop 
narrow cores. 



functions of the compactness of the neutron star and of 
its spin frequency. (The exact forms of the coefficients 
a,2n can be found in Morsink et al. 2007.) They find 
that even the shapes of the most rapidly rotating stars 
can be accurately described by truncating the series at 
N = 2. Since the Hartle-Thorne metric is formally of 
second order, we terminate the series in Equation (I7| at 
N = 1. 

In order to model the appearance of a neutron star 
to an observer at infinity, we use the ray-tracing code 
detailed in Baubock et al. (2012). In this code, an im- 
age plane is located at a large distance from the neutron 
star where the spacetime is asymptotically flat. Photon 
trajectories originate on the image plane and are traced 
backward through the curved spacetime using a fourth 
order Runge-Kutta algorithm. If the photons intersect 
the neutron star surface, the algorithm returns the 6 and 
<t> coordinates on the stellar surface, the energy of the 
photon, and its angle with respect to the surface nor- 
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Fig. 2. — Line profiles computed by taking different properties 
of spinning neutron stars into account. In each case, the star is 
rotating at 700 Hz. As in Figure [I] the stars have a radius of 
10 km, a mass of 1.4 Mq, and an inclination of 20° . The three lines 
show the results of calculations in three different approximations: 
the Schwarzschild+Doppler approximation for a spherical star, the 
Schwarzschild+Doppler approximation for an oblate star, and the 
Hartle-Thorne approximation. The last case is identical to the 
700 Hz profile shown in Figure [T] 
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Fig. 3. — Line profiles for four different inclinations to the ob- 
server's line of sight. The neutron star has a radius of 10km, a 
mass of 1.4 Mq, and a rotational frequency of 700 Hz. At high in- 
clinations, the profile displays only Doppler broadening. At lower 
inclinations however, a narrow peak appears. 



mal. Photons that pass beyond the star by ten percent 
of the distance to the observer are halted and discarded. 
For the purposes of this paper, we have assumed that the 
emissivity is constant over the surface of the neutron star 
and independent of the beaming angle of the photons. 

The ray-tracing code produces a map on the image 
plane of the ratio of the energy of the photon at its place 
of origin to its energy at infinity. In order to calculate 
the emission line profile for a distant observer, we use a 
Monte-Carlo algorithm to sample points across the image 
plane. At each point, we find the redshift by interpolat- 
ing the four nearest grid points of the ray-tracing simu- 
lation. Points falling in regions corresponding to photon 
paths that miss the neutron star are discarded. 

Additional effects at the neutron star surface may 
broaden the observed spectral lines beyond the purely 
kinematic and gravitational effects considered here. 
Mechanisms such as pressure broadening, thermal 
Doppler broadening, and Zeeman splitting are all ex- 
pected to contribute to the total width of observed lines. 
However, for moderately spinning stars such as those 
considered in this paper, the combined effect of these 
mechanisms is expected to be small compared to the ro- 
tational broadening (Ozel 2012). Therefore, we neglect 
such broadening mechanisms and consider hereafter the 
intrinsic widths of atomic features to be infinitesimal. 

3. LINE PROFILES 

Figure [l] shows the profiles of atomic features emitted 
from the surface of a neutron star with a radius of 10 km 
and a mass of 1.4 M Q , for different spin frequencies and 
observer inclinations. For each simulated profile, we nor- 
malize the line so that the bolometric luminosity is equal 
to unity. 

In the absence of rotation, the only effect the geometry 
of the star has on the line profile is to shift its observed 
energy by the gravitational redshift at the stellar surface. 
In the case of a 1.4 M Q star with a radius of 10 km, the 
energy of the line center is redshifted by a factor of 0.77. 

The zeroth order effect of rotation on the line profile 
is the broadening associated with the Doppler shift be- 
tween the approaching and receding edges of the star. 
At low spin frequencies, the width of the line is deter- 
mined almost entirely by this effect. A slight asymmetry 
appears in the profile that is due to the first order ef- 
fect of relativistic beaming — the photons emitted on the 
blueshifted portion of the star are beamed toward the 
observer, shifting the peak of the observed spectrum to 
the right as shown in Figure [T] The profiles correspond- 
ing to high observer inclinations closely match the results 
of, e.g., Ozel & Psaltis (2003), in which only first-order 
Doppler shift corrections were considered. 

At low inclinations we find a significant deviation from 
the line profile predicted by the Schwarzschild + Doppler 
broadening. In that simplified approach, the width of 
the line is determined by the line of sight velocity of 
the neutron star surface and scales with the sine of 
the inclination to the observer. We find, however, that 
neutron stars with relatively rapid rotations (such that 
their quadrupole moments are large) and low inclinations 
(such that the Doppler effects are small) generate line 
profiles with remarkably narrow cores that are strongly 
peaked toward the blue end of the spectrum. These peaks 
are evident in the higher frequency simulations shown in 
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Fig. 4. — Contours of constant redshift at the neutron star surface. The left panel shows the redshift for an oblate neutron star in the 
Schwarzschild+Doppler approximation, while the right panel includes our fiducial value for the neutron star as calculated by Laarakkers 
and Poisson (1999). In both cases, the star has a mass of 1.4 Mq, a radius of 10 km, an inclination of 20° to the observer's line of sight, 
and is rotating at 700 Hz. The star is rotating counterclockwise when viewed from above. In the right panel, the quadrupole moment of 
the neutron star causes a local maximum in the redshift to the right of the rotational pole. This maximum is responsible for the narrow 
core of the atomic features from moderately spinning neutron stars. 

the bottom panel of Figure [l] 

The deviation at low inclinations from the 
Schwarzschild+Doppler predictions is due primarily 
to the quadrupole moment of the neuron star, with 
a small additional effect from the oblate shape of the 
surface. To disentangle the relative roles of the Doppler 
shift, the quadrupole moment, and the oblateness, we 
simulated a neutron star with identical mass, radius, 
and spin frequency but made the surface spherical and 
performed the calculation in the Schwarzschild+Doppler 
approximation. The result is shown in the solid line in 
Figure [2] This profile shows only the expected Doppler 
broadening and relativistic beaming. When we allow for 
the shape of the star to become oblate, the profile shifts 
toward the red, as shown with the dashed line. This 
shift is caused by the increased gravitational redshift 
in the region near the pole. This effect is overpowered, 
however, by the contribution of a nonzero quadrupole 
moment. As shown in the dot-dashed line of Figure [2l 
the quadrupole causes a narrow peak in the line profile 
that overshadows the effect of oblateness. 

This narrow peak vanishes at higher inclinations, as 
shown in Figure [3] The peak is very strong at inclina- 
tions below 30° , but the effect of the extra quadrupole 
decreases rapidly as inclination increases because of the 
increasing relative importance of the Doppler broaden- 
ing. When viewed perpendicular to the axis of rotation, 
the peak vanishes and the profile can be described purely 
by Doppler broadening. 

The appearance of a narrow peak in our simulated 
spectra arises from the fact that the addition of a 
quadrupole moment causes a significant portion of the 
stellar surface to have roughly the same effective red- 
shift. This effect is a combination of the varying Doppler 
shift across the surface and the quadrupole compo- 



nent of the gravitational redshift. We show in Fig- 
ure E] contours of constant redshift on the stellar sur- 
face for two stars: the right panel depicts a neutron 
star with our fiducial value of the quadrupole moment, 
while the left panel shows the same calculation performed 
in the Schwarzschild+Doppler metric for an oblate neu- 
tron star. On the star with no quadrupole, the energy 
decreases monotonically toward the right (the receding 
side) due to the changing Doppler shift. In the other 
case, there is an additional effect due to the quadrupole 
that is symmetric about the pole of rotation. The result 
is a local maximum of energy that is offset slightly from 
the pole. The energy of this local maximum corresponds 
to the energy of the peak in Figure [l] 

In Figure [5j we present another view of the combi- 
nation of effects that give rise to these particular line 
profiles. The top panel shows the total redshift as a func- 
tion of colatitude for a constant value of the azimuth cf> 
close to the local maximum for the configuration shown 
in the right panel of Figure [4] For a spherical star with 
no quadrupole, the redshift changes monotonically across 
the surface due to the Doppler shift. An oblate star with 
no quadrupole moment is more redshifted in the vicin- 
ity of the pole, causing the shift seen in Figure [2] In 
contrast, for an oblate star with the fiducial quadrupole 
moment, a local maximum appears that is offset from 
the pole. 

By taking various ratios of the curves, we can identify 
the effect of the oblateness and the quadrupole moment, 
as shown in the lower panel of Figure [5] For comparison, 
we also show the quadrupole component of the gu term 
of the metric. The similarity between the redshift ratio 
and the metric curve strongly suggests that the redshift 
maximum is caused by the quadrupole moment. 

4. CONCLUSIONS 
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Fig. 5. — The top plot shows the redshift on the neutron star sur- 
face as a function of colatitude 9 for a fixed value of the azimuth 
<j>, for the neutron star parameters depicted in Figure 111 The three 
lines correspond to three different configurations of the neutron 
star: one in which the star is spherical and has no quadrupole mo- 
ment, one in which the star is oblate with no quadrupole moment, 
and one with both oblateness and a quadrupole moment. In the 
lower panel, the solid line shows the ratio of rcdshifts on a spheri- 
cal neutron star to those on an oblate one. The dotted line shows 
the ratio of the redshifts of a neutron star with a quadrupole mo- 
ment to those on a star without. The dash dotted line shows the 
quadrupole contribution to the gtt component of the metric. 

Second order effects of the rotation of neutron stars 
introduce significant corrections to the observed profiles 
of atomic features from their surfaces. The dependence 
of these profiles on the inclination of the observer devi- 
ates strongly from the expected sin(# ) behavior that is 
characteristic of Doppler effects. The principle cause of 
this deviation is the extra quadrupole mass moment in- 
duced by the rapid rotation of the star. The combination 
of the quadrupole moment and the Doppler shift due to 
the motion of the neutron star surface cause a strong 
and narrow feature to appear in the simulated spectra of 
emission lines. This narrow core is strongest at relatively 
low inclinations (6>o < 30°) and high spin frequencies 
(Jrts > 500 Hz). 

Figure [6] shows a contour plot of the simulated full 
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Fig. 6. — Contour plot of the full width at half maximum of the 
line profile from a spinning neutron star of 1.4 Mq and a radius of 
10 km as a function of spin frequency and inclination angle. 

width at half-maximum of emission lines for a range of 
spin frequencies and inclinations. At high inclinations, 
the interpretation of line widths is straightforward, with 
narrow lines corresponding to slow rotation and broad 
lines implying fast rotation. At lower inclinations, how- 
ever, there is often a degeneracy. Narrow lines could be 
caused by slow rotation or by the higher quadrupole in- 
duced by rapid rotation. Clearly, line broadening can 
only be used as a reliable measure of neutron star radii if 
additional information is available on the inclination of 
the observer and the spin frequency. 

The presence of a narrow core in atomic features also 
increases the range of inclinations at which these fea- 
tures are detectable. At high spin frequencies and rela- 
tively low observer inclinations, atomic features from re- 
alistic stars are narrower than those expected from purely 
Doppler effects (Chang et al. 2006). This may allow for 
such features to be discernible from continuum spectra 
and translates into a larger solid angle over which they 
are detectable 

The models presented here may shed new light on ob- 
servations of the source EXO 0748-676. XMM observa- 
tions of this source showed evidence for narrow absorp- 
tion lines (Cottam et al. 2003), which were later used 
to to constrain the mass and radius of the neutron star 
(e.g., Ozel 2006). However, Lin et al. (2010) argued that 
the subsequent detection of the 552 Hz spin frequency 
and the large amplitude of burst oscillations were incom- 
patible with the narrow observed width of the absorption 
lines and concluded that these lines did not originate at 
the neutron star surface. 

In their calculations, Lin et al. (2010) use the 
Schwarzschild+Doppler approach to model the line spec- 
tra. As shown above, however, the inclusion of second- 
order effects leads to narrow profiles for much higher in- 
clinations at a given spin frequency. Therefore, it appears 
plausible that the narrowness of the absorption lines is 
not incompatible with the high amplitudes of the ob- 
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served burst oscillations and, therefore does not exclude 
the possibility that they were emitted from the surface. 



A detailed analysis of the EXO 0748-676 case will be 
presented in a future paper. 
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